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Gauge invariant correlators in QCD are studied on the lattice. A systematic determination of the correlation lengths for gluon field 
strength correlators and quark correlators is made. The measurement of the gluon and quark condensates is discussed. 


1 Introduction 


The gauge invariant gluon field strength correlators are 
defined as 

V^.pA^.C) = (0|T(G;,(a.)5:^,-c(x,0)G^,(0)) |0) (1) 

with 5'“^^ q{x, 0) the parallel transport from 0 to x along 
the path G, in the adjoint representation 

Sadj.cix, 0) = Pexp (^i Apiy)Tadj (2) 

T^dj 8''^® the group generators in the adjoint representa¬ 
tion. 

Dp,vpa{x,C) depends on the choice of G: in what 
follows we will take for G a straight line, and drop the 
dependence on G. 

Higher correlators are usually defined by parallel 
transport to a fixed point xq. 

Fermion correlators are defined as 


^,(x) = (0|T {^{x)Sfund{x, O)MV(O)) |0) (3) 

M* is a generic element of the Clifford algebra of the 7 
matrices. Sfund is the analog of the transport (^) in the 
fundamental representation, and again we have once and 
for all assumed for the path G a straight lina. 

By use of general covariance argument£r0 Dpupu can 
be parametrized in terms of two independent invariant 
form factors V{x^) and Pi(x^) 


D. 


pypa(x^ — {gpipPi/a 9p.a-gup) ) V Dl{^X )j -f (4) 


pXpPi/fj X pX,y Pi/p XpXi/g^fj -\- Xi/XfjPpp^ 


dVi 

dx"^ 


Similarly one can prove that all the correlators i) vanish 
by T, P invariance, except the correlator Sq correspond¬ 
ing to M* = I, the identity matrix 

S'o(x) = (0|r [■ip{x)Sfundix,0)')pi0)) |0} (5) 

and the vector correlator (M® = 7^), with 7^ in the 
direction of x 

^5'v(x) = (0|T (v;(x)S'/w(a:,0)7'®V'(0)) |0) (6) 

The physical interest of the above correlators stems from 
the basic idea of the ITEP sum rulefl the long distance 
modes of QCD are described by a slowly varying back¬ 
ground, made e.g. of instantons, on which high momen¬ 
tum perturbative fluctuations are superimposed. In the 
O.P.E. low modes generate the condensates, and high 
frequency modes the corresponding coefficient functions 
G„(x) 

r(j;.(x)>(0)):^^G„(x)O„ (7) 

n 

= G/(x)/ + Cg{x)GpAQ)GpA0) + 

+ C,!, (x) ^ mf^f (0)V’/ (0) -b ... 

As is well known expressing the left hand side in terms 
of a dispersive integral, relates masses and widths of res¬ 
onances in e“'"e“ —*■ hadrons to the condensates 

G2 = ®(0|G;,(0)G“,(0)|0) , (0|m/^/(0)^/(0)|0) 

(SVZ sum rules). 
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By use of this idea it was proposed in refQid that 
the gluon condensate G 2 could be determined from the 
spectrum of bound states of heavy QQ systems. If the 
correlation length of the slow varying field, A, is much 
bigger than the typical time of the bound system, then 
its effect is in all respects a static Stark effect on the 
levels, and the gluon condensate can be extracted from 
it. 


A more detailed analysis involves 'D^^pa-{x), and 
gives a shift depending on the parameter 


p = X 


mqol 

4 


( 8 ) 


where 'ilmqol is the typical time of the low lying levels 
of the system, and A is the correlation length defined as 


V{x) ~ G '2 exp(—a;/A) (9) 

|tc|^oo 

Measuring A was the motivatioiii to investigate for the 
first time (x) on the latticfl. The computation was 

done in quenched SU(2) and gave a surprisingly small 
value of A 

Acs 0.16 fm (10) 

shaking the very bases of the SVZ approach. 

A stocl^stic model of the vacuum was subsequently 
developecu’El, in which observables are expressed in terms 
of invariant field strength correlators, and a cluster ex¬ 
pansion is made. The basic assumption of the model is 
that higher order clusters are negligible. 

The quark correlator S'o instead, known also as “non 
local fermion condensate” enters in the construction of 
the wave functions of hadrons, iii particular in the com¬ 
putation of the pion form factOB. 

A technical breakthrough, the use of cooling or 
smearing procedures to polish short distance fluctua¬ 
tions, leaving long distance physics unchangeca, allowed 
a better determination of correlatorfijO’Ej. 

The physical motivations to study correlators on the 
lattice are in conclusion 


(i) understanding A, and the basis of SVZ sum rules. 

(ii) measuring condensates from first principles. 

(iii) more generally providing inputs from first princi¬ 
ples to the community of stochastic vacuum prac¬ 
titioners. 


2 Cooling-smearing correlators 


in replacing a link by the sum of the inverse “staples” 
attached to it 




z 


◄ 


A 


( 11 ) 


Since the action density is 


► 


5 ~ ^ M j withn^i,= 

( 12 ) 

this procedure makes locally S = 0. In the euclidean 
region S plays the role of energy and the replacement (0 
locally minimizes S, whence the terminology “cooling”. 

Like any local procedure cooling rit times affects dis¬ 
tances d by a diffusion process, with 



(f nt 


(13) 


According to eq. (El), rit can be made sufficiently large 
to eliminate short range fluctuations but not enough to 
modify long range correlation^. Fluctuations are thus 
reduced by orders of magnitude without changing long 
distance physics. 

Gauge invariant correlators xan be represented on 
lattice by the following operators 


V 


L 

fiixpcr 



D □) 


A series expansion in a gives 

V^p,pM) cs Z^a^Vp,p.{d) + 0{a<^) (14) 

By the cooling procedure 0{a^) terms disappear and pos¬ 
sible renormalizations Z of the field Gpi, tend to 1. So 
finally 

Vli,p,{d) ~ a^Vp,p,{da) (15) 

In order to use cooling profitably the distance d has to 
be ~ 3 — 4 lattice spacings at least. At a given /3 this 
corresponds to some physical length Imin- If we want 
Imin to be small, say 0.1 fm, since the lattice must be 
at least 1 fm across, the lattice size must be (lOZmm)^- 
Going to small distances requires big lattices. 


3 Correlators, OPE and renormalons 

The SVZ sum rules are based on the OPE of the corre¬ 
lator 


Short range fluctuations in lattice configurations can be 
smoothed off by a local cooling procedure which consists 
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n;..(g) = J d^xe^<^^0\T{jpix)M0)) |0) (16) 

= n(g^) - q^qu) 










The OPE gives 


10 ' 


n(g^) ~ ciJ + C 2 ^ + . .. 


(17) 


The first term corresponds to the perturbative expan¬ 
sion. That expansion, however, is not Borel summable, 
and is ambiguous by terms of order with ^ the 

renormalization scale (Renormalons). As a consequence 
the second term in eq.(p^ is intrinsically undefined. This 
is a basic and unavoidable drawback of perturbation the¬ 
ory, reflecting the. fact that perturbative vacuum is not 
the ground statsid. 

However, keeping the first few terms in the per¬ 
turbative expemsion of ci and C 2 gives a consistent 
phenomenologjEl. 

The same happens for the correlators. The OPE of 
the invariant form factors has the form 


V{x^) + C 2 G 2 + (18) 


The second term is again undetermined by renormalons 
coming from the perturbative expansion of ci. As in the 
SVZ sum rules we shall assume that the first few terms of 
the perturbative expansion work, and use it to determine 
G 2 . We shall parametrize the lattice determination of 
'D{x^) and as 


^Vl{x^) = + Aoe-l"l/^ (19) 

\x\^ 

Lv^j^\x^) = ( 20 ) 


Eq.(|l^, (^^ obey the OPE eq.(^, and reflect the ex¬ 
istence of a mass gap in the theory. 


In ref0 the first term of eq. dH), ([!), ( [ 20 | ) was com¬ 
puted in perturbation theory and subtracted. The resid¬ 
ual term was an exponential and A could be extracted 
from it, giving A = 0.16 fm. 


In refM’B quenched SU{3) was studied. The above 
parametrization gave a good fit to the data (fig.l) with 
A = 0.22 fm and G 2 = (0.14 ± 0.02) GeV"*, a value 
larger by^n ordeiLof magnitude than the phenomenolog¬ 
ical valufi. In recj full QCD with 4 staggered fermions 
was studied, at quark ma^s arriq = 0 . 01 , arug — 0 . 02 . 
The results are in this cas£3 


io‘ 



Fig.l = V + 'Di+x^dVi/dx'^ and D±fa* = 

V -\-Vi versus x. The lines correspond to the best 
fit to eq.(19),(20). 


arUg = .01 A = (.34 ± .02) fm G 2 = .015 ± .003 GeV'‘ 
arUg = .02 A = (.29 ± .02) fm G 2 = .031 ± .005 GeV^ 


A common feature to all the determinations is that 
All ~ Ao/10. In full QCD the correlation length is big¬ 
ger, and could agree with the basic philosophy of reffl. 
Also the value of the condensate G 2 is smaller than the 
quenched value and agrees with phenomenology. 


By use of the 


relatioi 0 


d 

dirif 


G2 


24 - 
Oq 


60 = 11 - -Nf 


one can extrapolate in to/ to the physical value of G 2 
getting 


G 2 ~ 0.022 ± .006 GeV'^ ( 21 ) 


in agreement with sum 


rules 


determinatioiS. 


Similar arguments allow to extract G 2 from the mea¬ 
surement of the average value of the density of action 
(plaquette): a^in the level of rigour is the same as for 
SVZ sum ruletl, at least if only a few terms are kept 
of t^^ ^|rt^urbative expansion of the coefficients in the 


A detailed analysis of the behaviour of the correlators 
at finite temperature, in the idcinity of the deconfining 
transition Tc was made in refO. There 0(4) invariance 
is lost reducing to 0(3) and 5 independent form factors 
exist. 


The main result is that magnetic correlators are un¬ 
changed across Tc, while electric correlators have a sharp 
drop (fig.2,fig.3). 
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Fig.2 The electric longitudinal correlator versus 
distance, for different values oiTjT^. 



Fig.3 The magnetic longitudinal correlator versus 
distance, for different values of TjTf.. 

The scalar quark correlator 

5o(x) = -(0|T(^(x)5^^(0))|0) (22) 



to eq.(23). 

^=0.16 ±0.04 = 0.26 ±0.01 (24) 

A similar determination at amq = 0.02 /3 = 5.35 full 
QCD, /3 = 5.91 quenched, where a ~ 0.12 fm, give indis¬ 
tinguishable results: 

^=0.26 ±0.04 = 0.37 ±0.01 (25) 

Putting the two determinations together gives 

A/m^ = 1.5±0.3 (26) 

The typical correlation lengtlcis now ^ l/m^ in agree¬ 
ment with the approach of relij. 

A determination of Aq and a study of its relation to 
the quark condensate (ipip) is on the way. 

4 Discussion 


has recently teen determined in full QCD and in 
quenched QClia. 

A comparison has been made between the determi¬ 
nations in full QCD, at given values of the quark masses 
arug = 0.01, arug = 0.02 and of the lattice spacing a, and 
in quenched QCD at the same values of these physical 
parameters. No difference has been found, within errors, 
indicating that quark loops do not affect appreciably the 
quark correlator. 

A sensible parametrization for the lattice regulator 

Sq{x) = a^Ao exp(-x/A/) ± (23) 

Simulations have been performecS at arUg = 0.01 in full 
QCD at P = 5.35 and quenched QCD at /3 = 6.0, which 
correspond both to a lattice spacing a ~ 0.10 fm. No 
appreciable difference is found and in both cases (fig.4) 


The typical correlation length is rather small for gluon 
correlators: A = 0.16 fm for quenched SU{2), 0.22 fm for 
quenched SU{S), 0.32 fm for full QCD with 4 flavours. 
It is bigger for fermion correlators where 

A/ ~ (1.5±0.3)to"^ 

Condensates can be determined, despite the presence 
of renormalons, by the same philosophy used for SVZ 
sum rules. G 2 = (0.022 ± 0.006) GeV"' in full QCD is 
consistent with the determination by sum rules. It is 
an order of magnitude bigger in quenched QCD, G 2 = 
(0.14 ±0.02) GeV^. 

The behaviour of correlators at the deconffning tran¬ 
sition is consistent with expectations. 

Our determinations are useful to phenomenology and 
to test models of QCD vacuum. 
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